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We develop a theory of coherent quantum oscillations in two, in general interacting, qubits measured 
continuously by a mesoscopic detector with arbitrary non-linearity and discuss an example of SQUID 
magnetometer that can operate as such a detector. Calculated spectra of the detector output show 
that the detector non-linearity should lead to mixing of the oscillations of the two qubits. For 
non- interacting qubits oscillating with frequencies fii and 0.2, the mixing manifests itself as spectral 
peaks at the combination frequencies Q.i±Q,2. Additional nonlinearity introduced by the qubit-qubit 
interaction shifts all the frequencies. In particular, for identical qubits, the interaction splits coherent 
superposition of the single-qubit peaks &tQ,\ = 0,2. Quantum mechanics of the measurement imposes 
limitations on the height of the spectral peaks. 



I. INTRODUCTION 

Quantum measurements represent an important part 
of quantum information processing, quantum computing, 
and evolution of quantum systems in general. Counter- 
intuitive features of the measurement process related to 
the collapse of the wave- function of the measured system 
continue to attract interest to the "problem of quantum 
measurements" . Physics of mesoscopic solid-state qubits 
and detectors provides convenient tools for studying this 
problem in the most interesting case of quantum sys- 
tems that are large on atomic scale - see, e.g., chapters 
on quantum measurements in^ . Many features of quan- 
tum measurements manifest themselves directly in the 
regime of "continuous" measurements in which the mea- 
sured system evolves in time being continuously affected 
by the detector back-action which implements the wave- 
function collapse dynamically. Simple example of this 
regime is provided by the linear weak measurements of 
coherent quantum oscillations in one qubit^^^ which have 
been demonstrated experimentally in^. An interesting 
property of the continuous weak measurements is that 
the spectrum of the detector output has features that 
characterize directly the quantum mechanics of measure- 
ment and quantum nature of the qubit oscillations. 

One of the suggested tools in quantum-information ap- 
plications of quantum measurements are quadratic mea- 
surements, realized when the detector response to the in- 
put signal is quadratic. Quadratic measurements should 
make it possible to monitor products of operators of dif- 
ferent quantum systems, and can be used, for instance, to 
implement simple schemes of error-correction* , or to en- 
tangle non-interacting qubits^. The purpose of our work 
is to study the regime of continuous quadratic measure- 
ments concentrating on the case of two qubits which is 



the simplest system that reveals non-trivial character- 
istics of quadratic measurements. Currently, systems of 
two coupled mesoscopic qubits are almost routinely stud- 
ied in experiments - see^'^^""'^''. As discussed in this work, 
for two qubits, quadratic measurements are equivalent to 
measurements with arbitrary nonlinearity. Qualitatively, 
the main feature of continuous non-linear measurements 
is the mixing of coherent oscillations in individual qubits 
that for non-interacting qubits leads to appearance of 
spectral peaks at combination frequencies f2i±r22, where 
VLj are the frequencies of individual oscillations. In con- 
trast to mixing of classical oscillations, only these two 
combination frequencies appear in the detector output, 
and as in the linear regime, intensity of all spectral peaks 
is limited by the quantum mechanics of measurement. 

The paper is organized as follows. Section II intro- 
duces a model of mesoscopic solid-state detectors used in 
this work and gives an explicit description of one prac- 
tical detector which realizes this model and can operate 
in the purely quadratic or non-linear regime. Section 
III describes the two-qubit system and derives explicit 
equations for the evolution of the density matrix of this 
system in the measurement process. In Section IV we 
use these equations to calculate output spectral density 
of the non-linear detector measuring the two-qubit sys- 
tem in several situations. Section V provides concluding 
remarks. 



II. MESOSCOPIC QUANTUM MEASUREMENTS 

In this Section, we describe a generic model of quantum 
measurements with a mesoscopic solid-state detector and 
provide a detailed discussion of one specific example of 
such a detector. 
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A. General model of a mesoscopic detector 

Although mesoscopic detectors can have quite differ- 
ent physical implementations and include, e.g., quantum 
point contacts (QPC)^^"-^'''^^"^^, normal and supercon- 
ducting SET transistors^^"^^, SQUID magnetometers^^ 
and generic mesoscopic conductors'^^^'^^, the operating 
principle of all these detectors is essentially the same. 
Measured quantum system controls, through an opera- 
tor X, the transmission amplitude t{x) of particles (they 
can be electrons. Cooper pairs, or magnetic flux quanta) 
between the two reservoirs. The flux of these particles 
provides then the information on the state of this sys- 
tem. Schematics of the detector of this type measuring 
two qubits is shown in Fig. 1. 
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FIG. 1. Diagram of a mesoscopic detector measuring two 
qubits. The qubits modulate amplitude t of tunneling of de- 
tector particles between tlie two reservoirs. 

Quantitatively, the Hamiltonian of the detector-system 

coupling for such a detector consists of detector tunneling 
modulated by the measured system and can be written 
as 



(1) 



where ^, are the detector operators that create excita- 
tions when a particle is transferred, respectively, forward 
and backward between the reservoirs. For the quantum- 
point-contact (QPC) detector, which represents the sim- 
plest realization of the model discussed in this Section, 
^, describe excitation of electron-hole pairs in the QPC 
electrodes. 

We make several additional assumptions about the de- 
tector. We suppose that the tunneling between the de- 
tector reservoirs is weak and can be accounted for in the 
lowest non-vanishing order in the tunneling Hamiltonian 
(1). Under this assumption, the precise form of the in- 
ternal detector Hamiltonian is not important and dy- 
namics of measurement is deflned by the correlators: 
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dti^Hm, (2) 



where the angled brackets denote averaging over the de- 
tector reservoirs which are taken to be in a stationary 
state with the density matrix po- {■■■) = TtdI-'-Pd}- 
The correlators (2) set the scale r± = 2Re7± of the for- 
ward and backward detector tunneling rates. The corre- 
lators (C(t)C)) that do not conserve the number 



of tunneling particles are assumed to vanish. As can be 
seen more explicitly from the example of the SQUID de- 
tector discussed below, this condition holds even in the 
case of superconducting detectors, where such "anoma- 
lous" correlators can exists in general. Under the condi- 
tion of large bias voltage important in our model of the 
detector operation, the Cooper-pair tunneling is incoher- 
ent, and anomalous correlators indeed vanish. 

Another assumption is that the characteristic time of 
the detector tunneling is much shorter than that of the 
evolution of the measured system. One of the conditions 
implied by this assumptions is that the energy bias AE 
for tunneling through the detector, which sets one of the 
tunneling time scales, AE~^, is much larger than the typ- 
ical energies Eq of the measured system. In the example 
of the QPC detector, AE = eV, and this condition means 
that the bias voltage V across the QPC is sufficientlj^ 
large. For short tunneling times, the functions ^{t), ^^{t) 
in Eq. (2) are effectively ^-correlated on the time scale of 
the dynamics of the measured system. Condition of the 
large energy bias AE for the detector tunneling leading 
to the correlators Eq. (2) being ^-correlated is the natural 
part of the measurement model: it enables one to neglect 
quantum fluctuations in the detector in the frequency 
range that corresponds to that of the measured system, 
and makes the detector response in this frequency range 
classical. 

Combined with the assumption of weak tunneling, van- 
ishing correlation time in the correlators (2) makes it 
possible to write down simple evolution equations for the 
density matrix p of the measured system. Indeed, the 
time evolution of p in the interaction representation with 
respect to the tunneling Hamiltonian (1) is given by the 
standard expression: 

p{t)^TTD{SppDS^}, S = TeM-if dt'Hrit')}. 



(3) 

If the detector operators in Eq. (2) are (5-correlated, one 
can keep only the "non-crossing diagrams" in the pertur- 
bation expansion of Eq. (3) in Ht- Evolution of p{t) in 
Eq. (3) is governed then by the following equation: 



p = T+Ppt + rJpP - {'y+W + 7-Fi)p- 



(4) 



Equation (4) describes the measurement-induced part of 
the evolution of an arbitrary measured system within our 
generic model of a mesoscopic detector. 

One general remark that should be made here is that 
for some detectors, e.g. SET transistors, there are 
regimes of operation, when the particle transfer through 
the detector consists of more than one steps and can not 
be characterized by one transmission amplitude. Our 
measurement model is not applicable in these regimes. 
We focus on the case of one-step transfer, however, since 
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only in this case the detector can be quantum-hniited: 
intermediate steps introduce additional back-action de- 
phasing for the measured system without increasing in- 
formation contained in the detector tunneling rate. 

Evolution (4) of the density matrix p of the measured 
system is reflected in the detector output: the particle 
current between the detector reservoirs. Using the same 
logic that lead to Eq. (4), one can see that the assumed 
largo difference between the time scale of detector tun- 
neling and evolution of the measured system makes it 
possible to reduce the expression for the particle current 
to the form of the operator I in the space of the measured 
system: 

I = {T+-T_)tH. (5) 

This equation gives both the dc current 

(/) = Tr{7po} , (6) 

where po is the stationary solution of Eq. (4), and the 
current spectral density 

/>oo 

Si = So + 2 (ircoswr(Tr{/e^"[/po]} -(/)'). (7) 
Jo 

Here Sq represents the usual noise (typically, a mixture 
of shot and thermal noise) associated with tunneling: 

So = {r+ + r_)TT{thpo}, (8) 

and e^"^ [A] denotes the evolution of the matrix A during 
time interval r governed by Eq. (4). 

B. DC SQUID as the non-linear detector 

Before discussing applications of our measurement 
model that is expressed quantitatively by Eqs. (4)-(8), 
we give one explicit example of the mcsoscopic detec- 
tor (different from the simplest QPC detector that was 
already mentioned several times in this Section) that re- 
alizes this model. The detector is the strongly-biased dc 
SQUID shown schematically in Fig. 2. It consists of two 
Josephson junctions with critical currents 7i,2 that are 
included in small superconducting loop shunted by a re- 
sistor with some impedance Z{ui), assumed to be feature- 
less, Z{lo) = Z{0) = R, at low frequencies of variations 
of the measured magnetic flux through the SQUID 
loop. When the loop inductance is small, the difference 
between the two Josephson phases yi,2 across the two 
junctions is directly linked to $: 

Vi-V2 = 27r$/$o = 4> , 

where $o = '"'^i/e is the magnetic flux quantum. In this 
regime, the SQUID is equivalent to a single Josephson 
junction, with the effective amplitude of Cooper-pair tun- 
neling controlled by the flux Due to interference of the 



Cooper-pair timneling in the two SQUID junctions, the 
total tunneling amplitude is equal to the sum of tunnel- 
ing amplitudes in the two Josephson junctions, with the 
amplitude in each junction being proportional to its crit- 
ical current Ij . The total amplitude can then be written 
as 

where ii = [Ixjl'i)^!'^ and = characterize the 
asymmetry of the two junctions. 

The operating principle of the SQUID as a detector 
coincides with the one discussed in the beginning of this 
Section: variations of the flux <p around some bias point (j) 
lead to variations of the Cooper-pair tunneling amplitude 
(9) affecting the transfer rate / of Cooper pairs through 
the SQUID. Cooper-pair tunneling rate is reflected in the 
detector output: deviations V of the voltage across the 
SQUID, V = -2eIR, from the value Vo = Rio induced 
by the dc current bias Iq (see Fig. 2). An important 
feature of the SQUID detector is that in the case of iden- 
tical junctions: ii = ^2 = 1, at the bias points when 
^/27r is integer, the tunneling amplitude t{(f>) (9) varies 
quadratically as a function of (f>, with vanishing coeffi- 
cient of the linear response. As shown explicitly in the 
next Section, under this bias condition the SQUID can 
act as the purely quadratic detector. 
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FIG. 2. DC SQUID detector formed by two Josephson 
junctions with critical currents 7i,2 in a small superconduct- 
ing loop shunted by the impedance Z{u)). The SQUID is bi- 
ased by the dc current 7o, the voltage V is the measurement 
output and reflects the variations in the rate of incoherent 
Cooper-pair transfer through the SQUID with amplitude t 
controlled by the flux <I>. Existence of the two interfering con- 
tributions to the tunneling amplitudes makes it possible for 
the SQUID to operate as a quadratic detector. 

The Cooper-pair tunneling Hamiltonian in the SQUID 
can be expressed through the tunneling amplitude (9) 
similarly to Eq. (1): 

(10) 

where (fi{t) = {fiit) + (f2{t))/2 is the average Joseph- 
son phase across the SQUID which includes fluctuating 
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component accumulated due to the equilibrium voltage 
fluctuations across the impedance Z{ll!). The fluctua- 
tions couple the Cooper-pair tunneling to Z{uj) so that 
the transfer of a Cooper pair across the SQUID creates 
electromagnetic excitations in this impedance, the pro- 
cess described by operators e^*"^^*^ in Eq. (10). In this 
way, impedance Z{u}) provides dissipation necessary to 
make the Cooper-pair transfer irreversible. 

When the dc bias current Iq and associated voltage 
Vo = RIq across the dc SQUID are sufficiently large, one 
can treat the tunnel Hamiltonian (10) as perturbation 
(i.e., the dc current through the SQUID is much smaller 
than the bias current Iq) and the system dynamics at low 
frequencies can be described as incoherent tunneling of 
Cooper pairs through the SQUID^'^. This regime satisfies 
both of the conditions on the detector dynamics discussed 
in the general detector model, and when the control flux 
(j) is created by the measured quantum system, measure- 
ment dynamics is governed by the same general equation 
for the system density matrix (4), where the correlators 
(2) are now given by the following expressions'^^: 
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(11) 



Here Rq = nh/Ae'^ is the "quantum resistance". Un- 
der the conditions of weak linear coupling of the SQUID 
detector to the measured systems, the correlators (11) 
that give the rates of incoherent Cooper-pair tunneling 
are the only relevant properties of the SQUID operating 
as a linear quantum detector'*. If the variations of the 
external flux </> controlling Cooper-pair tunneling are not 
small, the SQUID can act as non-linear, and in partic- 
ular, purely quadratic detector. Dynamics of the mea- 
sured system is reflected in the Cooper-pair current in 
the SQUID in exactly the same way as in the general 
model [see Eqs. (5)-(8)], and is converted into the detec- 
tor output voltage V as discussed above. 



III. TWO QUBITS MEASURED CONTINUOUSLY 
BY THE NON-LINEAR DETECTOR 

The main focus of this work is on the measurement of 
the system of two, in general interacting, qubits. The 
qubit Hamiltonian is: 



(12) 



i=i,2 



where Aj is the tunnel amplitude and ej is the bias of 
the j-th qubit {j = 1,2), v is the qubit-qubit interaction 
energy, and cj's here and below denote the Pauli matri- 
ces. The qubits are coupled to one detector (see Fig. 1) 
through their basis- forming variables a^, and the con- 
trol operator x modulating the detector tunneling (1) is 



X = cial+ C2(T^. It gives the following expression for the 
tunneling amplitude t: 



i{x) = to -h ^ Sjal + Xcrlal 



(13) 



The last term in this expression appears due to non- 
linearity of the dependence of the transmission amplitude 
t on variable x. If the linear terms in Eq. (13) vanish, 
6j — 0, while A 7^ 0, one has purely quadratic detector. 
In the case of two qubits, since the Pauli matrices sat- 
isfy the condition = 1, Eq. (13) represents the most 
general dependence of t on cr|.^^ In the example of the 
SQUID detector discussed in Sec. II, when the two flux 
qubits are coupled to it, the normalized flux </> through 
the detector is: 



+ 2 



where Scpj characterizes the strength of coupling to the 
7th qubit; and the average flux 4> sets the detector op- 
erating point. With such a coupling, the amplitude (9) 
of the Cooper-pair tunneling in the detector is given by 
Eq. (13) with 

to = [iie*^/^ -I- 126-'^^^] cos5^!>i cosS(t>2, , 
6j = i[he'^l^ - i2e-^'^/2] cos 5<j)f sin Scpj , f ^ j , 
A = -[iie''^/2 ^ i^e-'^/"^] sin sm6(j)2 • 

If = 2?™ with integer n, and the SQUID is symmetric: 
ii = i2 = l,thc linear coupling coefficients Sj vanish and 
the SQUID has only quadratic response. 

For qubit-detector coupling of the form (13), it is con- 
venient to write the general equation (4) for the density 
matrix p of the two-qubit system in the "measurement" 
basis of eigenstates of the cr^ operators, | "ft) , | "f J,) , 1 11) , 
and I ii). Each state \k) of this basis is characterized by 
the magnitude tk of the transmission amplitude (13): 

ti — to + 61 + 62 + X , t2 — to + 61 ~ 62 — X , 
ts = to — ^1 + ^2 — X, ti = to — Si — 62 + X , 

and associated value of the detector tunneling current, 
for which we use the obvious notations: 





= (r+ 


-r_)|to + <5i + (52 + A|2 


hi 


= (r+ 


-T_)\to + 6i-S2-Xf 


hi 


= (r+ 


-T_)\to-5i + S2-X\^ 


hi 


= (r+ 


-r_)\to-di-62 + Xf 



Combining measurement-induced evolution (4) with the 
evolution due to the qubit Hamiltonian Ho we get the 
equation for p in the measurement basis: 



Pki = -IkiPki - i[Ho,p]ki , 

7fei = (i/2)(r+ + r_)|ife-tH'. 



(14) 
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In Eq. (14), the Hamiltonian Hq includes two renormal- 
ization terms: 

Ho^ Ho + SH + SH' . 

The first one is due to imaginary parts of the correlators 

(2): 

SH = Y,Sejai + 6ualal (15) 
j 

where 

5ej = Re((5jto + <^j'A*)lm(7_ + 7+) , 
= Rc((5i(52 + toA*)Im(7_ + 7+) , 

and = 1, 2, with j' ^ j. The second term 6H' is due 
to phases ipki = avg{tkt*) of the transfer amplitudes tk 
and is defined by the following relation: 

[6H',p\ki = (r+ - r_)\tkti\ sinipkipki ■ (16) 

Note that 5H' can not always be cast in the form (15) in 
which it can be absorbed in the renormaHzation of qubit 
energies in the Hamiltonian (12). This can be done if 
the difference between the tunneling amplitudes tk are 
small, \tk — ti\ <C \tk\- In this work, we will assume that 
this condition on 6H' is satisfied, and the renormalized 
Hamiltonian has the same form (12). 

Equation (14) describes the time evolution of the 
qubits averaged over different measurement outcomes, 
which in our measurement model arc represented by the 
number n of particles tunneled through the detector. Be- 
cause of this averaging, the qualitative effect of measure- 
ment in Eq. (14) is the "back-action" dcphasing of differ- 
ent states of the measurement basis with the rates jki ■ In 
general, Eq. (4) contains information about the dynamics 
of n and can be used to write down the evolution equa- 
tions conditioned on the specific measurement outcomes 
(see Ref. 9). The averaged equation (14) is sufficient for 
calculation of the output spectrum of the detector which 
is the purpose of this work. 

The assumption of short detector tunneling time that 
lead to Eq. (14) makes this equation valid even in the 
regime when the detector-qubit coupling is strong and 
the dephasing rates "/ki arc large in comparison with the 
rates of evolution due to the Hamiltonian Hq. Straight- 
forward numerical solution of Eq. (14) combined with 
Eqs. (7) and (8) gives the spectrum of the detector out- 
put. (All numerical plots of the spectra are obtained be- 
low in this way.) An example of such a spectrum for most 
generic set of parameters of the qubit-detector system is 
shown in Fig. 3. For these parameters all six intervals 
between four energy levels of the two-qubit Hamiltonian 
(12) are different, and show up as six finite-frequency 
peaks in the spectrum of the detector output. There is 
also the zero-frequency peak that reflects the dynamics 
of detector-induced transitions between the energy lev- 
els. Qualitatively, the width of all peaks corresponds to 



the rates of transitions between the energy levels, with 
finite frequency peaks are also broadened by the "pure 
dephasing" terms in the evolution equation for the den- 
sity matrix of the system [see Eq. (19) below] that are 
not related to the transitions. 

Although so far the nonlinear mesoscopic detectors of 
the type considered in this work have not been realized 
experimentally, many elements of our detector model 
were indeed demonstrated. For instance, the QPC de- 
tector, which is the simplest realization of this model, 
was used to measure coherent oscillations in a quantum- 
dot qubit^*. Parameters of the detector-qubit system 
used in Fig. 3 and other numerical results presented below 
(in particular the characteristic detector tunneling rate 
r+|foP on the order of the qubit oscillation frequency A) 
are consistent with those in experiments. 
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FIG. 3. Output spectrum of a nonlinear detector measur- 
ing two qubits with "the most general" set of parameters. 
Six peaks in the spectrum at finite frequencies correspond to 
six different energy intervals in the energy spectrum of the 
two-qubit system. The zero-frequency peak reflects dynamics 
of transitions Ijctwccn energy levels. Detector parameters are: 
Si = 0.1, S2 = 0.07, A = 0.09 (all normalized to to). In this 
Figure, and in all numerical plots below we take r+|to|^ = Ai, 
r_ — 0, and assume that the detector tunneling amplitudes 
are real. 

In the situation without any symmetries in the qubit 
parameters, matrix elements of the detector-qubit cou- 
pling are non-vanishing between all four energy eigen- 
states of the two-qubit system (12). Moreover, the as- 
sumption of the short detector tunneling time implies 
that the transition rate is independent of the energies of 
these states: i.e., if one views the detector as the reser- 
voir producing dephasing for the measured system, the 
effective reservoir temperature is much larger than the 
system energies. This means that the stationary density 
matrix of the system is: 

Po = //4, (17) 

where / is the unity matrix in the four-dimensional space 
of the two-qubit system. This means that the dc detector 
current (6) in this regime is: 
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(r+-r_)(|io|' + |-5i|^ + |52|' + |A|^), 



(18) 



and associated noise So (8) is given by the same expres- 
sion with r+ — r_ replaced by r+ + r_. If, however, 
the qubit parameters have some symmetry (examples are 
given in the Sections that follow), matrix elements sat- 
isfy "selection rules" and vanish for transitions into some 
states. The stationary density matrix po is then propor- 
tional to the unity matrix confined to a smaller subspace 
of the Hilbert space of the two-qubit system. In this case 
the dc current (/) and the noise 5*0 are given by equations 
different from Eq. (18), and the number of peaks in the 
detector spectrum is less than the maximum of 6 peaks. 

It is interesting to note that the selection rules that 
determine the number of finite-frequency peaks in the 
output spectrum of the detector allow sometimes for a 
simple qualitative classical interpretation. For instance, 
as will be demonstrated more rigorously below, in agree- 
ment with the case of classical oscillations, if the qubit- 
qubit interaction vanishes, and the detector is purely lin- 
ear, the output spectrum contains at most two peaks 
that correspond to oscillations in the individual qubits 
with frequencies r2i,2- If the detector is non-linear, the 
spectrum acquires two more peaks at the combination 
frequencies fli ± $72- The maximum total number of 
peaks is also equal to four if the detector is linear but 
the qubit-qubit interaction is finite. 

Quantitatively, the detector spectrum can be obtained 
from the solution of Eq. (14) for evolution of the qubit 
density matrix. Despite the relative simplicity of its 
numerical solution, Eq. (14) can be solved analytically 
only in a few cases [see, e.g., Eq. (24) below]. More de- 
tailed analytical results can be obtained if the dcphas- 
ing rates jki are small in comparison to the intervals 
Wnm = En — Em between eigenenergies of the Hamil- 
tonian Hq (12). In this limit, it is convenient to trans- 
form Eq. (14) into the basis of eigenstates \n) of Hq, 
Haln) = En\n), where one can separate components of 
the density matrix p that evolve with different frequen- 
cies. Then, making use of the standard approach (see, 
e.g.,'^"'') equation for p can be simplified by neglecting the 
terms that mix these components. Written in the energy 
eigenstate basis, Eq. (14) is reduced in this way to the 
following form: 

Pnm = —itJ^nmPnm + {^+ + T-) X 

p^m P^n 
^nm / ^ Ppp I ^mp 

P + (1 - Snm) X] Ppq^^i^ipkm) ■ (19) 

p {p,q) 

where tpq are the matrix elements of the tunneling am- 
plitude t (as operator in the qubit space) in the basis of 
energy eigenstates, and the last sum is taken over the 
pairs {p, q) of states that satisfy the "resonance" condi- 
tion: Ep — Eq = En — Em , {p,q) 7^ {n,m), and there- 



fore contribute to the same spectral peak of the detector 
output. Solution of Eq. (19) can be used to calculate the 
average detector current and its spectrum using the same 
Eqs. (5)-(8). In the next Section, we use both Eq. (14) 
and Eq. (19) to describe the time evolution of the two- 
qubit system and resulting detector output in different 
situations. 



IV. SPECTRAL DENSITY OF THE DETECTOR 
OUTPUT 

Dynamics of the detector-qubit system and spectral 
density of the detector output depend on several char- 
acteristics of the system, most important of which arc: 
the degree of non-linearity of the detector-qubit cou- 
pling, symmetry of parameters of individual qubits, and 
strength of the qubit-qubit interaction. In this Section, 
we calculate the spectra of the detector output in several 
regimes that differ in terms of these characteristics. Sub- 
section A discusses the case of purely quadratic coupling 
of the detector to qubits that, in general, are different and 
interacting. Subsection B deals with the detector-qubit 
coupling which includes both the linear and quadratic 
terms but only for identical qubits. In the last Subsec- 
tion C, we consider the most general case of different 
and interacting qubits measured by the detector with ar- 
bitrary non-linearity limiting ourselves to the regime of 
weak detector-qubit coupling. 



A. Purely quadratic detector 

We begin the discussion of the detector output spec- 
trum by considering the purely quadratic measurement 
of unbiased qubits (Sj = ej =0). In this case, dynamics 
of measurement governed by Eq. (14) is such that the 
Hilbert space of the two-qubit system is split into two 
two-dimensional subspaces D± with no transitions be- 
tween them, so that the system evolves independently in 
each subspace. If the basis of states in the subspace D± 
is chosen as 

D± ^ {^(1 TT) ± I ii)) . it) ± I Ti))} > (20) 

the Hamiltonian (12) is split into two independent parts 
H± acting within D±: 



H± = -(l/2)[(Ai ± A2)a, + va,] . 



(21) 



(By putting the bars on the Pauli matices a in this eciua- 
tion, we want to distinguish these Pauli matrices act- 
ing in the subspaces D± from those describing individual 
qubits.) 

Qualitatively, the qubit dynamics within each subspace 
can be viewed as oscillations between parallel and anti- 
parallel configurations. The oscillation frequencies are 
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equal to Ai ± A2 for non-interacting qubits and are 
shifted upwards by interaction which creates energy dif- 
ference between the parallel and anti-parallel configura- 
tions. Since the purely quadratic detector distinguishes 
these two types of qubit configurations, the current op- 
erator (5) within each subspace in the basis (20) is: 



/=(/)+ laa. 



(22) 



where the average detector current (/) is given by 
Eq. (18) with 5j = 0, and 

h = (r+ - r_)2Re(toA*) = /tT - ^Ti 

has the meaning of the amplitude of the detector current 
oscillations reflecting the qubit oscillations between the 
parallel and anti-parallel configurations. 

Quantitatively, the evolution equation (14) within D± 

is: 



0, 







(23) 



where F = 2(r+ -I- r_)|Ap is the measurement-induced 
dephasing of the basis states (20) of each subspace. For 
non-interacting qubits, solving this equation and calcu- 
lating the spectral density of the detector output accord- 
ing to Eq. (7) we get^: 



2/2(Ai ± A2)2F 



(a;2 - (Ai ± A2)2)2 + F2a;2 ' 



(24) 



where 5o = (F+ + F_)(|toP + |A|2). 

The two spectral densities (24) correspond to two pos- 
sible outcomes of measurement: the qubits found in one 
or the other subspace D±, the probability of the out- 
comes being determined by the initial state of the qubits. 
Each of the spectral densities coincides with the spectral 
density of the linear detector measuring coherent oscilla- 
tions in one qubit^. Similarly to that case, the maximum 
of the ratio of the oscillation peak versus noise 5*0 for 
each spectrum Sf{u) is 4. As one can see from Eq. (24), 
this maximum is reached when the measurement is weak: 
|A| ^ |to|, and the detector is "ideal": arg(toA*)=0; and 
only F_|_ or F_ is non-vanishing. If, however, there is 
small but finite transition rate between the two subspaces 
that mixes the two outcomes of measurement, the peak 
height is reduced by averaging over the two spectral den- 
sities (24). This situation is illustrated in Fig. 4 which 
shows the output spectra of the purely quadratic detec- 
tor, when the subspaces D± are mixed by small qubit bias 
e. Since the stationary density matrix (17) is equally dis- 
tributed over all qubit states, the two peaks of the spec- 
tral densities (24) are mixed with equal probabilities, and 
the maximum of the ratio of the oscillation peak heights 
versus noise So for the combined spectrum Si{lo) is 2. 
Spectrum shown in Fig. 4 for £ = 0.1 Ai (solid line) is 
close to this limit. 
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FIG. 4. Output spectra of a purely quadratic detector 
measuring two non-interacting qubits. Small qubit bieis 
El = £2 = e (solid line) creates transitions that lead to averag- 
ing of the two main peaks at combination frequencies Ai ± A2 
[see Eq. (24)]. Further increase of e (dashed line) makes ad- 
ditional spectral peaks associated with these transitions more 
pronounced. The strength of quadratic qubit-detector cou- 
pling is taken to be A = O.lSto- 

For interacting qubits, analytical results for the spec- 

trmn can be obtained in the case of weak back-action 
dephasing F <C |Ai ± A2I, when the evolution equations 
(23) within each subspace can be transformed into the en- 
ergy eigenstates representation and simplified there sim- 
ilarly to Eq. (19). In this limit, the shape of all spec- 
tral peaks is Lorentzian, with the decay rate 7 of the 
off-diagonal elements of the density matrices p^^^ in the 
eigenstate representation, 

7=(F++F_)|A|2(l + ^.Vf]2^), 

determines the width of the peaks at the combination 
frequencies shifted by interaction: 



n± = [(Ai±A2)2 



21 1/2 



(25) 



In the case of non-vanishing interaction strength v, the 
eigenstates of the Hamiltonians (21) have different am- 
plitudes of states (20) with parallel and anti-parallel 
qubit configurations and therefore have different values 
of the average detector current (22) . This means that for 
r <C |Ai ± A2I, nonvanishing interaction v creates the 
spectral peak at zero frequency which reflects the transi- 
tions between the energy eigenstates. The width of this 
peak 7o is determined by the rate of transitions between 
these states: 

70 = 2(r+ -t- F_)|A|2(Ai ± A2)VJ1^ 

The total spectra Sf{uj) within each subspace D± are: 



Sf{u)=So 



10 



'-(w - n±)2 -1-72 



-)]) • 



(26) 
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Finite qubit bias should lead to averaging of the two spec- 
tra (26) similar to that discussed in the case of non- 
interacting qubits and illustrated in Fig. 4. 

If the difference of the qubit tunneling amplitudes be- 
come small, |Ai — A2I <^ F, qubit dynamics in the D_ 
subspacc can no longer be viewed as coherent oscillations 
between the basis states (20), but rather as incoherent 
transitions between these states. The rate of these inco- 
herent transition t^^ can be foimd by treating Ai — A2 
as perturbation in the evolution equation (23) for the 
density matrix p^~^ : 



l (Ai 
2 



A2)2r 



+ r2 



(27) 



Since the basis states (20) are characterized by different 

values (/)+ Jq of the detector current - sec Eq. (22), these 
transition give rise to the spectral peak in the detector 
output spectrum at zero-frequency: 



57(w) = 5o + 



1 + T2a;2 



(28) 



Equation (28) for the zero-frequency peak in the detector 
output is valid for arbitrary relation between the back- 
action dephasing rate and interaction strength u. It re- 
produces both the zero-frequency peak in the spectrum 
(26) when T <^ v, and if F » u, the peak in ob- 
tained from Eq. (24) in the limit |Ai - A2I < F. If 
the qubits arc identical, Ai = A2, the transitions (27) 
are suppressed and both basis states of the subspace £>_ 
represent separate outcomes of quadratic measurement. 



B. Non-linear detector measuring identical qubits 

In the case when the two qubits have the same set of 
parameters, it is convenient to discuss dynamics of the 
two-qubit system using the language of the total "spin" 
S = (cTi +a2)/2. The qubit Hamiltonian (12) and the de- 
tector tunnel amplitude (13) providing the qubit-detector 
coupling can be written in terms of S like this: 



Hf) — —sSz 
t = to + 26Sz 



A5, 
-X{2Sl-l) 



(29) 
(30) 



where e, A, 5 without indices denote the same quantities 
as for individual qubits. The state (| |t) — I Ti))/V2 with 
S = Q does not evolve in time under the Hamiltonian 
(29) and represents one of the measurement outcomes 
characterized by the dc detector current 

(j) = (r+-r_)|io-A|2 

and flat output spcctrmn Si(uj) = (F+ + r_)|fo — Ap. 

Three other, S = I, states arc mixed by measurement 
and represent the second measurement outcome. Insert- 
ing Eq. (30) into Eq. (5) and using the fact that in the 



S = 1 subspace = Sz, we get the following expression 
for the detector current operator in this subspace: 



/ = (r+ - F_)|to - A|2 + aiSz + '2a2S\ 



2 

Z 1 



(31) 



where ai,2 have the meaning of the amplitudes of current 
oscillations between different qubit states: 



ai 
02 



4(r+ 

2(F+ 



.)(Re[ioA1 + W) = (/tt 



^ii)/2, (32) 
+ /xi-2%)/4. 



Similarly to Eq. (17), the stationary density matrix 
within the 5 = 1 subspace is po = 1/3- Taking the aver- 
age over the three eigenvalues of Sz operator, Sz — 0, ±1, 
with this density matrix, we see that the dc detector cur- 
rent in the measurement outcome that corresponds to the 
5 = 1 subspace is 



F_ 



■[2(|to|' + |A|2) + |io + A|2 + 8|(5|% (33) 



{!) = 



and can be written as (/) = {I^ + + /||)/3. 

To calculate the spectral density of the detector current 
(31), we consider first non-interacting qubits and limit 
ourselves to the case of weak measurement which can 
be described conveniently by going to the basis of the 
energy eigenstates. Using the vector intuition for spin 
operators, one sees directly that an appropriate rotation 
of the Hamiltonian (29) brings it to the form 



H 



-nSz, o = (A^-F£^) 



2x1/2 



(34) 



with three energies {— Q, 0,0}. Upon this rotation, the 
operator Sz in the tunneling amplitude (30) changes ac- 
cordingly 



Sz^isSz + AS,)/ft. 



(35) 



In the regime of weak measurements, spectral peaks 
of the detector output at different frequencies are de- 
termined by the evolution of different groups of matrix 
elements of the operator s = pol, which evolve indepen- 
dently one from another. A peak at some finite frequency 
Lu is determined by the off-diagonal matrix elements of s 
between the states with the energy difference equal to w. 
This means that the total number of the matrix elements 
relevant for a given uj coincides with the number of times 
this interval occurs in the energy spectrum. In the situ- 
ation considered in this subsection, s = I/'S, where the 
current operator I is given by Eq. (31) in which Sz is 
transformed according to Eq. (35). The structure of the 
energy levels of the Hamiltonian (34) implies that there 
should be three peaks in the spectrum of the detector 
output: at w :^ 0, 2f2, and the zero-frequency "relax- 
ation" peak. 

If the basis states are numbered in the direction of 
increasing energy, H\l) = H\2) = 0, etc., the 

peak at a; ~ 217 is determined by the matrix element 
Si3 which satisfy the simple equation that follows from 
Eq. (19): 
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si3 = {i2n - 7)si3 , (36) 

and the initial condition Si3(0) = /is/S, where 

7i3 = asAVn' 

is the matrix element of the cmrent operator given by 
Eqs. (31) and (35). The decoherence rate 7 in Eq. (36) 
is: 

7=1 [T12 + T23 + 2Tis + (r+ + T-){4s\6\/nf] , 

where the last term represents the pure dephasing, and 
Tij is the rate of transitions between the states i and j: 



T12 = 2(r+ + T_)\6n + Xel^A^/n" , 

^l3 = (^+ + ^_)|A|2AV^^^ 



(37) 



and the rate is given by the same expression as 
with e — > — £. All the rates are obtained from Eq. (19) 
in which the matrix elements of the operator of the the 
tunnel amplitude are given by Eqs. (30) and (35). Equa- 
tion (36) means that the spectral peak at a; ~ 20 has 
simple Lorentzian form: 



1 



3 (w - 20)2 + ^2 



(38) 



and the background noise coincides with Eq. (33) for 
the dc current in which r_|_ — r_ is replaced by r+ + r_. 

The structure of the a; ~ H peak of the spectral density 
of the detector output is determined by the time evolu- 
tion of the matrix elements S12 and S23 which satisfy the 
coupled system of equations following from Eq. (19): 



S12 
S23 



iOsi2 

iOs23 



^lSl2 + KS23, 
C2S23 + KS12. 



(39) 



Initial conditions in these equations are the same as in 
Eq. (36), Sij(O) = /ij/3, where the current matrix ele- 
ments are: 

7i2 = ai(A/\/20) + a2(\/2eA/02) , 

and 723 is given by the same expression with £—>—£. 
Again, the relaxation rates ^„i, m = 1,2, and the rate k 
of the "transfer of coherence" are obtained by combining 
the matrix elements of the transmission amplitude (30), 
(35) with Eq. (19): 



U = {T+ + T-)[\5\\2 
x4Re(<5A*)(e/0)3 + |Ap(An 

K = (^+ + ^_)(|5po2- 



^ aVo2) - (-1)™ 

e^A^ + 2e4)/l]4] ^ 
-|A|V)(2AV0^). 



Solving Eq. (39) by diagonalization of the matrix of evo- 
lution coefficients, we see that the spectral peak at w ~ O 
consists in general of two overlapping Lorentzians with 
different line- widths and amplitudes: 



m=l,2 



(a; -0)2 + 72 



(40) 



where 



7m = (6+6)/2 + (-l)™-D, 

[(6-6)V4+«^]^/^ 

r2 , r2 (-1)™ 



D 

A-rr 



(41) 



D 



[2hhK-{ll 



II 



where for later convenience we introduced notations: 
li = Ii2, h = l23- In the case of unbiased qubits, e = 0, 
the situation becomes much simpler: = ^2, -'^12 = -^23, 
and one of the the Lorentzians vanishes, A2 = 0. In this 
case the peak at a; ~ has the form of one Lorentzian 
with the line- width 71 = (r+-|-r_)(|(5|2-|-|A|2) and ampli- 
tude Ai = 2a\. For purely linear measurement: A = 0, 
\5\ <^ \to\, this implies that the maximum peak height 
2af/37i is hmited by (32/3)50.3*^ 

The spectral peak at w ~ is determined by the transi- 
tions between the energy eigenstates with the rates (37). 
The average current Ijj in the j-th state is different from 
the dc current (33) through the detector: 



7ii - (7) = ai(£/0) + a2(2£2 _ A^)/3n'^ , 
I22 - (7) = 2a2(A2 - 2£2)/302 , 



(42) 



733 - (7) = -ai(£/0) + a2(2£2 _ A^)/m'' 



Because of nonvanishing differences Ijj — (!) , transitions 
between the eigenstates generate low-frequency noise in 
the detector current. Quantitatively, solving the stan- 
dard kinetic equation for the evolution of occupation 
probabilities of the three eigenstates due to transitions 
(37) and substituting the solution into Eq. (7) we see 
that similarly to the peak at ~ O, the zero-frequency 
peak of the spectral density of the detector output con- 
sists of two overlapping Lorentzians with different widths 
rim and amplitudes Bm- 



Sj{u;) = So + l E J 

rn=l,2 



iVr, 



(43) 



Parameters of the two Lorentians are determined by the 
transition rates (37) and the currents (42): 



= ri2 
[r?2 



r23 + ri3-(-i)'"7?. 



r^3 + r 



13 



ri2ri3 — ri3r23 

a?(e/0)2 + Cm'aliA^ - 2eY/^^^ 



ri2r23 - 

,2/ A 2 



1/2 



1ln 

D 

7?m = |CmaT(e/iM" + C;,, 
+(-l)'"(ri2 - r23)2aia2£(A2 - 2e^)/n''] /D , 

Cm = 2D + (-i)'"(ri2 + r23 - 2ri3) , 



where m' is defined as m' = 1, 2, m' 7^ m. 

The situation again simplifies drastically for unbiased 
qubits, £ ^ 0, when only one of the Lorentzians has 
non-vanishing amplitude B = 20.2/3. The width of this 
non-vanishing Lorentzian is 



9 



T] = 3(ri2 + r23)/2 = 6(r+ + T_)\sf . 

Equations (38), (40), and (43) describe completely the 
output spectrum of the detector weakly coupled to iden- 
tical non-interacting qubits. In agreement with classical 
intuition, the spectral peak at twice the frequency CI of 
the oscillations in individual qubits appears only if the 
detector-qubit coupling is effectively non-linear, when ei- 
ther A 7^ or is not too small. The peak at individual 
qubit frequency Q has generically larger amplitude than 
the peak at 2U and can be viewed as the result of coherent 
superposition of oscillations in two qubits. (Note that the 
maximum peak height, (32/3)5*0, at a; ~ is larger that 
two times the height, 45*0, of the spectral peak associated 
with the oscillations in individual qubits.) This interpre- 
tation is supported by the behavior of the two peaks as a 
function of weak qubit- qubit interaction v which we con- 
sider now. Including the interaction term in the Hamilto- 
nian (29) in the evolution equation for the density matrix 
in the basis of energy eigenstates (34) one can see that 
the interaction has no effect on the dynamics at ~ 251 
until interaction strength v becomes comparable to 17. 

In contrast to this, the dynamics at w ~ is affected 
by much weaker interaction on the order of the detector- 
induced dephasing. The interaction breaks coherence be- 
tween oscillations in the two qubits and eventually splits 
the spectral peak at a; ~ O in two as the interaction 
strength v becomes larger than dephasing. We describe 
this splitting quantitatively limiting ourselves to the case 
of unbiased qubits, when O = A. In the regime of weak 
interaction ^ A, we can use the "rotating-wave" ap- 
proximation by keeping only those interaction terms in 
the evolution equations that do not mix components os- 
cillating with different large frequencies on the order of 
A. In this case, the peak at w ^ A is governed by the 
dynamics of the same matrix elements Si2, S23; as for 
non-interacting qubits, and Eqs. (39) for the dynamics 
of these elements now are: 
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FIG. 5. Evolution of the output spectrum of the non-linear 
detector measuring two identical unbiased qubits with the 
strength u of the qubit-qubit interaction. The qubit-detector 
coupling constants 5i,2 are taken to be slightly different to 
average the spectrum over all qubit states. The three solid 
curves correspond to j^/A = 0.0, 0.1, 0.2. In agreement with 
Eqs. (45) - (47), the peak at a; ~ A is at first suppressed and 
then split in two by increasing v, while the peak at u; ~ 2A is 
not changed noticeably by such a weak interaction. Dashed 
and dotted hnes show the regime of relatively strong interac- 
tion: iz/A = 0.5 and i^/A = 1.0, respectively, that is described 
by Eqs. (49) and (50). 

Equations (44) have the same structure as Eqs. (39), 
and for weak interaction, v/2 < k = 2(r+ + r_)|i5|^, 
result in the same form of the peak of the spectral density 
at w ~ A as in Eq. (40): two overlapping Lorentzians 

m = 



with different line-width 7„ 
1,2: 



and amplitudes A„ 



7rr 



A,fn = fl]^ (1 - 



V4] 



1/2 



(45) 



If u/2 K, then D 0, and the two Lorentzians com- 
bine to give the peak of the following form: 



(44) 



si2 = i(A -I- iy/2)si2 - Csi2 + KS23, 

S23 = i{A - U/2)S23 - ^S23 + KS12, 

where we took into account that for vanishing bias, 

Ci = 6 = (r+ + r.)(3|5p + |An^^ 



We see that for <C A, the effect of the interaction is 
just the shift of energy of the zero-energy eigenstate of 
the Hamiltonian (34) by 1^/2 relative to the two other 
eigenstates. 



5,(.)=50 + f (1-4) (^,^^^2+^2 ■ (46) 

For stronger interaction, i//2 > k, the two Lorentzians 
acquire different frequencies A ± D, D = [i^^/4 — k^]^^"^, 
while having the same line- width ^, so that the total spec- 
tral density at w ~ A is: 



i + K±{A-u)n/D 



(47) 



For v ^ K, spectral density (47) consists of two inde- 
pendent Lorentzians separated by the large frequency in- 
terval v. The amplitude of both peaks is a^/3^ and is 
limited by (16/9)5*0, the value that should be reached by 
an ideal hnear detector. 

Equations (45), (46), and (47) are valid in the regime 
of weak interaction, i/ ~ k <C A, when interaction affects 
strongly only the spectral peak at a; ~ A, which is first 
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broadened and then split in two with increasing interac- 
tion strength. To describe the detector output spectrum 
for stronger interaction ~ A 3> we need to diagonal- 
ize the Hamiltonian (29) with non-vanishing v. This can 
be done easily in the case of unbiased qubits e = Q, when 
the eigenstatcs of the Hamiltonian (29) coincide with the 
two eigenstatcs of the Hamiltonian (21) and the state 
(I TT) — I ii))/V2 of the subspace Z)_ introduced before 
in Eq. (20). The eigenenergies of these states are: 

{-f2+/2 , Z//2 , f2+/2 }, f2+ = (4A2 + i/^) V2 _ (43) 

The three finite-frequency peaks in the spectrum cor- 
respond to three energy intervals, ± v)/2 in the 
spectrum of these eigenstates. Since all energy intervals 
are different and therefore there is no transfer of coher- 
ence between them that exists, e.g., in Eqs. (39), the 
peaks are simple Lorentzians. Transforming the tunnel- 
ing amplitude (30) into the basis of eigenstates (48), we 
find all the rates in the evolution Eq. (19) and the ma- 
trix elements of the current operator (31) in this basis 
and find the parameters of these Lorentzians: amplitude 
A and the line- width 7, both defined as in Eq. (40): 

a; ~ n+ , A = 8(a2A/fi+)2 , 

7 = (r+ + r_)[2|<5|2 + (1 + ^Vf^+)|An , (49) 

Cc;~(fi+±i^)/2, A = a\{l^v/Vl+) , 
7 = (r+ + r_)[3(l T + (1 ± z^/0+)|A|2] . (50) 

One can see that parameters of the two lower-frequency 
peaks (50) agree in the regime A » 1/ k with the 
peaks in the spectral density (47), while the peak (49) at 
the largest frequency 0+ coincides with that described 
by Eq. (38) the case of unbiased non-interacting qubits. 

Figure 5 illustrates evolution of the output spectrum of 
the non-linear detector measuring identical qubits due to 
changing interaction strength. We see that this evolution 
agrees with the analytical description developed above. 
Weak qubit-qubit interaction ~ k ^ A suppresses and 
subsequently splits the spectral peak at w ~ while 
not changing the peak lo ~ 2fi. Stronger qubit-qubit 
interaction v c^i ^ ^ n shifts the uj c^i 20, peak to higher 
frequencies while moving the two peaks around a; ~ il 
further apart. 



C. Non-linear detector measuring different qubits 

As the last example of the output spectrum of the 
non-linear detector measuring two qubits we consider the 
most general situation when both the tunneling ampli- 
tudes and the detector-qubit coupling constants are dif- 
ferent for the two qubits. We assume that the qubit 
parameters are such that all energy intervals in the spec- 
trum of eigenstates are larger than the back-action de- 
phasing rate, i.e., the detector-qubit coupling is weak. 
We begin with the case of unbiased qubits, ej = 0, when 



the coherent oscillations in the two qubits should have 
the largest amplitude. At this bias point, the Hamilto- 
nian (12) breaks into two subspaces D± (20) and can 
be diagonalized directly. It is convenient to order the 
eigenstates taking into account these subspaces: 



{n+/2, -n+/2, n_/2,-n_/2} . 



(51) 



The energies 0± are defined in Eq. (25). The wavefunc- 
tions of the eigenstates numbered in this order are: 

IV'i) = [{a + m TT) + I ii)) + (/3 - a){\ U) + I iT))]/2 , 
\^2) = [{a - m TT) + I ID) + {a + m\ U) + I iT»]/2 , 
IV'a) = [{a + m TT) - I ID) + (/3 - m ID - I Ti))]/2 , 
IV'4) = [{a - m TT) - I ii)) + (a + ^)(| TT) - I Ti))]/2 , 
where 

a,/3=(l/y2)[l±(Ai+A2)/0+]i/2, 
aj= (l/\/2)[l ± (Ai - A2)/n_]i/2 _ 

Assuming that there are no "accidental" degeneracies, 
there are four different energy intervals in the spectrum 
(51) which should be reflected as four finite-frequency 
peaks in the detector output spectrum. For different 
qubit parameters, there are no selection rules and the 
detector back-action mixes all four state of the qubit 
system. The stationary density matrix is then given by 
Eq. (17) and the background detector noise So corre- 
sponds to the dc current (18). The largest (ri+) and 
the smallest (r2„) energy intervals occur only once in the 
spectrum (51), so that the peaks at w ~ 0± has the 
shape of simple Lorentians: 



UJ ~ fl± , Sjioj) = So 



1 



A, 



4{u;-n±)^+ii 



(52) 



Calculating the current matrix elements and all the 
rates in Eq. (19) for evolution of the density matrix in 
the basis of eigenstates (51) similarly to what is done for 
the description of decoherence in coupled qubits^^, we 
find parameters of the Lorentzians in Eq. (52). It is con- 
venient to express the peak amplitudes A± through the 
characteristic amplitudes of current modulation by the 
qubits analogous to the amplitudes ai,2 (32) for identical 
qubits. In the case of different strength of the detector- 
qubit coupling for the two qubits, we have three such 
amplitudes: 

an = {In + In - In - in)/"^ 

= 2(r+-r_)ReMi* + A5^], 

ai2 = (/TT+/iT-Ai-^ii)/4 (53) 
= 2(r+-r_)Re[to'52+A(5t], 

a2 = {In + In - In - in)/"^ 

= 2{T+-T_)Re[tQX*+6iS;]. 

For identical qubits, an and ai2 reduce to ai/2, while 
the definition of 02 in Eqs. (53) and (32) coincide. Qual- 
itatively, aij corresponds to the amplitude of modula- 
tion of the detector current due to oscillations in the j-th 
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qubit, and 02 - is the similar amplitude due to "collec- 
tive" oscillations of the two qubits between parallel and 
anti-parallel configurations. In terms of the amplitudes 
(53), A± are given by simple expressions: 

A± = 2a2[(Ai±A2)/f2±]2 (54) 

The line-widths 7± in Eq. (52) are: 

7± = (r+ + r_)[|5i|2 + \62\' + |A|2(i + i^'/nl)] . (55) 

In contrast to the energy intervals 0± which occur only 
once in the energy spectrum (51), the intervals (f2+ ± 
f2_)/2 occur twice each. This means that the spectral 
peaks at w ~ (0+ ± f7_)/2 are not simple Lorentzians, 
and their shape is controlled by the two matrix elements 
of the density matrix evolving according to a system of 
coupled equations identical with Eqs. (44). This means 
that these peaks consist of two overlapping Lorentzians 
each, and the output spectrum in their vicinity is given 
by the equation similar to Eq. (40): 

4„tr,2[c^-(0+±fi_)/2]2 + [7W]2 

for u> ~ (r2+ ± ri_)/2. Here the amplitudes A^^^ and 
line-widths are given by the same Eqs. (41), where 
now 

= (r+ + r_) [(1 ± ^^)|Ap + |<5ip + |^2p 



^(r+ + r_)[''^ 



^2 



{\Si\' + \62f) (56) 



±l)(|(5ip-|52p)] 



7W = (±1)-+^^[1t'''"^' + ^^^'/' 
V2 



+(T1)™^[1T 



Ainl/2 



In the case of non-interacting qubits and no quadratic 
coupling, u = X = 0, the peaks at (Jl+±fi_)/2 represent 
oscillations in the individual qubits: the peak at (Q_)_ — 
f2_)/2 describes oscillations in a qubit with a smaller A, 
while the peak at (f2+ -|- n_)/2 - in the qubit with larger 
A. Only one Lorentzian (to = 1) is then non- vanishing 
for each peak, and has the same parameters as it woiild 
have in the absence of the other qubit. For instance, if 
Ai > A2, we get from Eqs. (41) and (56) {ov u = X = 0: 



(r!+-r!_)/2 = A2 



(r+- 



r_)|5i 



2l^ A 



4+)=4a?„ 



This behavior is the natural consequence of the lin- 
earity of the detector-qubit system for z/ = A = 0. 



Non- vanishing qubit-qubit interaction and/or non-linear 
detector response modify parameters of these "linear" 
peaks according to Eqs. (41) and (56). In contrast to 
these peaks, the spectral peaks at Q± reflect directly 
the detector non-linearity and vanish in the linear regime 
when A = and 6j <C to- As one can see from Eq. (53), 
the characteristic amplitude 02 of these peaks scales for 
A = as a higher power of Sj than the amplitudes aij of 
the linear peaks. 

An example of the output spectrum of the non-linear 
detector measuring unbiased qubits with different tunnel- 
ing amplitudes is shown in Fig. 6. One can see that when 
the linear and non-linear coefhcient of the detector-qubit 
coupling are roughly similar, the linear peaks are more 
pronounced than the peaks at combination frequencies. 
Qubit-qubit interaction shifts all but the lower-frequency 
linear peak up in frequency and reduces both the am- 
plitudes of the higher-frequency peaks and the distance 
between them. 
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FIG. 6. Output spectra of the nori-lincar detector mea- 
suring two different unbiased (jubits. Solid line is the spec- 
trum in the case of non-interacting qubits. The two larger 
peaks are the "linear" peaks that correspond to the os- 
cillations in the individual qubits, while smaller peaks are 
non-linear peaks at the combination frequencies. Dashed line 
is the spectrum for interacting qubits. Interaction shifts the 
lower-frequency liner peak down and all other peaks up in 
frequency. Parameters of the detector-qubit coupling are: 
<5i = 0.12to, 52 = 0.09to, A = 0.08to- 

Analytical results of this subsection can be easily ex- 
tended to the finite qubit bias if the qubits are non- 
interacting. The eigenstates of the Hamiltonian (12) are 
then the products of the eigenstates of the individual 
qubits, so that the energy spectrum of the two-qubit sys- 
tem is: 

(i/2){i]i + ^2, ni - ^2, ^2 - ni, -{ni + ^2)}, (57) 

where ±nj/2 with D,j = (e? -|- A'jy^^ are the eigenener- 
gies of the j-th qubit. 

Similarly to the spectrum (51), there are four different 
energy intervals in (57). Two of them, fii ± ^2 = oc- 
cur only once and are reflected as simple Lorentzian peaks 
in the detector output spectrum. These "non-linear" 
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peaks represent the mixture of the individual qubit oscil- 
lations with frequencies f2i^2 and have non-vanishing am- 
plitude only if the detector response is non-linear. Quan- 
titatively, the output spectrum around this peaks can be 
written as in Eq. (52), where now the peak parameters 
are: 



A. 



AiA 



2\2 



7± = (^^ 



£2 I 2 
f2^' . 



(58) 



They are found by calculating both the rates in Eq. (19) 
for the density matrix, and the matrix elements of the 
current operator J: 



in the basis of states (57). In this expression, (/) is the 
average current (18). 

The energy intervals Qj occur twice each in the spec- 
trum (57). This means that the spectral peaks at a; ~ flj 
that correspond to the individual qubit oscillations have 
the same form as in Eq. (40), so that 



Si{uj) = So + 



E 

TO=1,2 



A. 



(^_0.)2+^^ 



(59) 



for a; ~ fij, where the amplitudes Aj^ and line- widths 
7jTO are given by Eqs. (41) with 

= (r+ + r_)[(i + [£,• )|5,- - (-i)™A(£,-v%')|' 

+(|A|2 + |5,.|')(A,-'M')'], 

= (r+ + r_)(A,Vf2,0'[l'5,f + \)^?{e,l^,f] , 
Ijm = [aij - {-ira2{ej,/nj,)] (A,/fi,-) , 

where j' ^ j. 

We see from Eq. (59) that for non-interacting qubits, 
the non-vanishing qubit bias just shifts the frequency po- 
sition of the liner peaks (56) without qualitatively chang- 
ing their shape. If both the bias and the qubit-qubit in- 
teraction are finite, the bias splits each of the linear peaks 
in two simple Lorentzians bringing the total number of 
the finite-frequency peaks in the spectrum of the detec- 
tor output to six as it should be in the generic situation 
(see, e.g.. Fig. 3). 



V. CONCLUSION 

To summarize, we have developed a theory of con- 
tinuous quantum measurements of coherent oscillations 
in two coupled qubits by a non-linear detector. Calcu- 
lated spectra of the detector output show that the detec- 
tor non-linearity leads to the appearance of the spectral 



peaks at the combination frequencies of the qubit oscilla- 
tions in the detector output. The spectra have the non- 
trivial dependence on the strength of the qubit-qubit in- 
teraction. For identical non-interacting qubits, the spec- 
tral peaks at frequency of individual qubit oscillations are 
superimposed coherently, with weak interaction break- 
ing this coherent superposition and splitting the oscil- 
lation peak. In general, qubit-qubit interaction should 
manifest itself qualitatively through the total number of 
peaks in the output spectrum: the total number of peaks 
is at most 4 in the case of non-interacting qubits, while 
it can reach 6 in the most general situation with both 
non-vanishing qubit-qubit interaction and the detector 
non-linearity. 
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